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We describe a method for determining the radiative decay properties of a molecule by studying 
the saturation of laser-induced fluorescence and the associated power broadening of spectral lines. 
The fluorescence saturates because the molecules decay to states that are not resonant with the 
laser. The amplitudes and widths of two hyperfine components of a spectral line are measured over 
a range of laser intensities and the results compared to a model of the laser-molecule interaction. 
Using this method we measure the lifetime of the A(v' = 0) state of CaF to be r = 19.2 ± 0.7 ns, 
and the Franck-Condon factor for the transition to the X(t) = 0) state to be Z = 0.987 ^'019 ■ 
In addition, our analysis provides a measure of the hyperfine interval in the lowest-lying state of 
A(v' = 0), A e =4.8 ± 1.1 MHz. 



PACS numbers: 33.50.Dq, 33.70.Ca, 33.15.Pw 

I. INTRODUCTION 

Atomic and molecular beams are often detected by 
laser- induced fluorescence (LIE). The technique offers 
high efficiency, excellent spectral resolution, and, if the 
beam is pulsed, high temporal resolution [ij. In pass- 
ing through a probe laser whose frequency is resonant 
with an electronic transition, the molecules scatter pho- 
tons from the laser into a detector. If the overall detec- 
tion efficiency and the mean number of photons scattered 
per molecule are known, the absolute molecular flux can 
be determined. The mean number of photons scattered 
per molecule is a function of laser intensity, but satu- 
rates when this intensity is high enough. In this paper, 
we study this saturation of laser-induced fluorescence in 
detail, and we use our results to measure the natural 
linewidth and the Franck-Condon factor of the A-X(O-O) 
transition in CaF. 

We distinguish two mechanisms for the saturation of 
laser-induced fluorescence. In the first mechanism, sat- 
uration occurs once the excitation rate R exceeds the 
spontaneous decay rate T. Further increasing the inten- 
sity does not increase the fluorescence, since this becomes 
limited by the rate of spontaneous emission. This first 
mechanism is most important for 'closed' systems where 
the excited state can only decay back to the initial state. 
In an 'open' system, which is more common, the excited 
state can decay either to the initial state or to one or 
more other states that are not coupled to the laser. In 
this case, the average number of photons per molecule 
cannot exceed 1/(1 — p), where p is the probability of 
returning to the initial state. The first saturation mech- 
anism still applies to this case, but now the fluorescence 
will also saturate when RT ^> 1/(1 —p), for then, in the 
time T taken to pass through the probe laser, the average 
number of scattered photons per molecule has reached its 
maximum value. Further increasing R by increasing the 



laser intensity cannot then increase the fluorescence. The 
first mechanism is treated in most text books on atom- 
laser interactions, e.g. 0, but it is the second mechanism 
that is responsible for saturating laser-induced fluores- 
cence in practically all atomic and molecular beam ex- 
periments that use cw lasers. The detection of molecular 
beams by laser-induced fluorescence is treated by Hefter 
and Bergmann Q, and in the earlier review articles by 
Kinsey Q and by Altkorn and Zare Q. 

When fluorescence saturates by either mechanism, the 
LIF spectral lines will broaden because the degree of sat- 
uration will be greater on the line centre than in the 
wings. This power broadening is very well known in the 
context of the first saturation mechanism, though again it 
is the second mechanism that causes power broadening in 
most high resolution LIF spectroscopy with atomic and 
molecular beams. By measuring the power broadened 
linewidth as a function of probe intensity the natural 
linewidth of the transition, and hence the excited state 
lifetime, can be determined, provided that other broaden- 
ing mechanisms are well controlled. It is also possible to 
determine the value of p by measuring the saturation of 
fluorescence with probe intensity. For an electronic tran- 
sition in a diatomic molecule, p will depend on rotational 
factors that are usually straightforward to calculate, and 
on the Franck-Condon factor for the transition which is 
often unknown. Thus, by measuring how the widths of 
the spectral lines broaden and how the amplitudes of the 
lines saturate as the probe intensity is increased we can 
determine lifetimes and Franck-Condon factors. These 
are important molecular parameters that are often rather 
difficult to measure with high precision. Their determi- 
nation has received renewed interest recently because of 
the possibility of laser cooling certain molecules where 
the lifetime is short and the Franck-Condon factor is very 
close to unity @. 

Lifetimes are usually measured directly by observing 
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FIG. 1: (Color online) Experiment setup 

the rate of fluorescence decay following excitation by a 
short laser pulse. For CaF, the lifetime of the A(v' — 0) 
state was measured to be 21.9 ± 4 ns by this method 0- 
Franck-Condon factors are sometimes measured directly 
by measuring the relative intensities of a sequence of vi- 
brational transitions, but more often are determined the- 
oretically using potential energy curves constructed from 
measured spectroscopic parameters. In the case of CaF, 
both parameters have been calculated recently using the 
multireference configuration interaction method with a 
large basis set [|[ , so our measurements provide a precise 
test of this advanced molecular theory. The technique we 
use here is a general one that could be applied to many 
other molecules. 



II. EXPERIMENT 

The experimental method is rather simple. The setup 
is similar to the one described in Q, and is shown in 
Fig.[T] A pulsed solenoid valve filled with a mixture of 
Ar (98%) and SF6 (2%) to a pressure of 4 bar, emits gas 
pulses of 100-200 /is duration at a rate of 10 Hz into a 
vacuum chamber maintained below 10~ 4 mbar. Immedi- 
ately outside the nozzle of this valve, at the moment when 
the gas density is highest, a Ca target is ablated by light 
propagating along x and produced by a Nd:YAG laser 
that provides pulses of duration 10 ns, energy 25 mJ and 
wavelength 1064 ran. This laser pulse defines the zero 
of time in the experiment. The Ca target is made in 
an inert environment by cutting strips of Ca, ~ 2 mm 
thick and ~ 4 mm wide, and then gluing them to the 
rim of a stainless steel disk. This disk is mounted in 
the vacuum chamber with its axis parallel to z, and with 
the Ca surface displaced a few mm from the beamline in 
the x direction. The flux of CaF decays, so after a few 
minutes of ablating a single spot of the target, it is ro- 
tated a little so that a fresh piece of Ca is exposed to the 
laser. The beam passes through a 2 mm diameter skim- 
mer, situated 70 mm downstream, into a second vacuum 
chamber maintained below 10 -7 mbar, and then through 
a detector situated 810 mm from the source. 



In the detector, the time-of-flight profile of the CaF 
molecular pulse is recorded, with a resolution of 10 /is, 
by laser-induced fluorescence at 606.3 nm. Light from 
the interaction region is imaged, with a magnification of 
2, onto the 28 mm diameter photocathode of a photo- 
multiplier tube operated in current mode. A rectangu- 
lar mask, 9 mm along x and 20 mm along z is placed in 
front of the photocatode to spatially filter the fluores- 
cence. This limits the visible part of the molecular beam 
to a 4.5 mm region along x, thereby reducing the Dopplcr 
contribution to the measured linewidth to about 5 MHz. 
This is important for an accurate determination of the 
natural linewidth. The probe laser beam comes from a 
tunable dye laser, passes first through an acousto-optic 
modulator (AOM), then through beam-shaping optics, 
then polarizing optics and finally propagates through the 
vacuum chamber (along x) via Brewster-angled windows. 
The AOM is used to modulate or scan the amplitude 
or the frequency of the light, the beam-shaping optics 
are used to produce a collimated circular beam approxi- 
mately 2 mm in diameter with an approximately top-hat 
intensity distribution, and the polarization optics ensure 
that the laser light is plane polarized along z. The power 
of the probe laser is monitored just before it enters the 
vacuum chamber, and the beam profile monitored by 
splitting off a portion of the beam onto a profiler whose 
distance from the beamsplitter is the same as that of the 
laser-molecule interaction region. 

The time-of-flight profile is used to determine the 
molecular speed distribution whose mean is typically 
600 m/s. Integrating the time of flight profile gives the 
total fluorescence, which we measure as a function of 
laser frequency to obtain a spectrum. We analyze laser- 
induced fluorescence spectra that contain just two spec- 
tral lines, the F = 1 and F = hyperfine components of 
the A-X(0-0)Qu(l/2) transition, separated by 123 MHz. 
A portion of the laser beam passes into an evacuated 
confocal cavity whose free spectral range and finesse are 
250 MHz and 2.1 respectively. The transmission of this 
cavity is recorded at the same time as the spectra, and 
this signal is used to linearize the spectral data and cali- 
brate the frequency axis. 

Figure [2] shows example spectra obtained at two dif- 
ferent probe intensities, along with typical cavity data. 
The line through the cavity data is a fit to an Airy func- 
tion, and the line through each spectrum is a sum of 
two Lorentzians, with the amplitude, centre and width of 
each Lorentzian floating in the fit. Because the Doppler 
broadening is very small in the experiment, this is the ex- 
pected lineshape for sufficiently low laser intensities. At 
higher laser intensities the spectral lines are power broad- 
ened and the exact lineshape can only be determined by 
the detailed modelling discussed in Sec. IIII1 Neverthe- 
less, the two-Lorentzian lineshape was found empirically 
to be a good fit to all the data used in our analysis. 
At very high laser intensities the lineshape departs from 
the two-Lorentzian model. In Fig. [21 the amplitude ra- 
tio between the F = 1 and F — lines is not 3, as it 
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FIG. 2: (Color online) Points, upper: typical laser-induced 
fluorescence spectra showing the two hyperfine components 
of the CaF A-X(0— 0) Qn(|) transition. The intensity at the 
centre of the laser beam was 11.5 W/m 2 for spectrum (i) and 
97.6 W/m 2 for (ii). Lines, upper: double Lorentzian fits to 
the data. Points, lower: transmission of a 250 MHz confocal 
cavity used to linearize the scan and calibrate the frequency 
axis. Line, lower: Airy function fit to the cavity data. 



would be from the degeneracy factor alone, but close to 
4.5 for both spectra. This happens because an excited 
F = 1 molecule is more likely to decay back to the orig- 
inal hyperfine state and scatter a second photon than is 
the case for an F = molecule. For spectrum (i), for 
which Iq = 11.5 W/m 2 , the full width at half maximum 
(FWHM) is 17 MHz for the F = line and 13 MHz for the 
F = 1 line. For the Iq — 97.6 W/m 2 data, the lines are 
strongly power broadened and the widths have increased 
to 39 MHz and 28 MHz. The power broadening is a direct 
result of the saturation in the scattered photon number, 
which in turn is due to the decay of excited molecules 
to states that are not resonant with the laser. Because 
F = 1 molecules can scatter more laser photons, the on- 
set of saturation occurs at higher intensity than for F = 
molecules, and so at a given intensity the F = 1 linewidth 
is smaller. 

Spectra such as those shown in Figure [2] were analyzed 
over a wide range of laser intensities. For each spec- 
trum, the scan was linearized, and then the amplitudes 
and widths of the two lines determined by fitting a two- 
Lorentzian model to the data. Fitting to Voigt functions 
instead of Lorentzians made a negligible difference to the 
fit residuals. 

For the majority of the data used in our analysis, the 
AOM was used to modulate the probe power from high 
to low on successive shots from the source so that we 
obtained two (nearly) simultaneous spectra, one at high 
power and the other at a lower power. The high value 
was fixed for the entire set of spectral scans, while the 
lower value was varied between scans. This allowed us 
to normalize all the scans so as to eliminate slow drifts 
in the molecular flux. In other data sets we used the 
AOM to modulate the frequency of the probe laser by 
123 MHz while scanning the frequency of the dye laser 



as usual. This gave us two (nearly) simultaneous spec- 
tra, one shifted by 123 MHz relative to the other, so that 
the F = 1 to F = amplitude ratio could be measured 
with high precision. Finally, some additional data were 
obtained by locking the frequency of the dye laser and 
scanning the probe frequency using the AOM. This has 
the advantage that the frequency scan is exactly linear, 
without the need for later correction using the cavity 
data. Since the bandwidth of the AOM is limited, only 
one hyperfine component could be scanned this way, and 
then only for low laser intensities where the linewidth is 
small. For these data, the rf power delivered to the AOM 
was controlled by a feedback loop so as to keep the laser 
power constant over the scan. 

III. MODEL 

Figure [3] shows the relevant molecular energy levels for 
the experiment. In the X 2 £ + (v — 0) ground state, the 
angular momentum of the rotating nuclei, TV, couples to 
the electron spin, S, to form a total electronic angular 
momentum J — N + S. Since S — 1/2, each rotational 
level (apart from N = 0) is split by the spin-rotation 
interaction into two components with J = N ± 1/2. Cou- 
pling to the 1—1/2 spin of the fluorine nucleus leads 
to further splitting of each level into hyperfine compo- 
nents labelled by F, the quantum number of the total 
angular momentum F = I + J . In the experiment, the 
laser drives Qn(l/2) transitions from the ground state, 
X 2 S+ (v = 0,N = 0, J=l/2), to the electronically excited 
state, labelled A 2 n 1/2 (v' = 0, J' — 1/2, /). Here, the or- 
bital angular momentum, L, is strongly coupled to the 
internuclear axis, and so too is the electron spin due to a 
strong spin-orbit coupling. The sum of their projections, 
fl = A + E, has magnitude 1/2, this being the spin-orbit 
component of lower energy. In the A state, the projec- 
tion of L onto the internuclear axis can either be A' = +1 
or A' = — 1. The energy eigenstates are the symmetric 
and antisymmetric superpositions of these two projec- 
tions. The degeneracy between these two eigenstates, 
which have opposite parity, is lifted by a small degree 
of mixing with other electronic states (A-doubling). The 
laser drives transitions only to the state labelled /, the 
one of negative parity, because the ground state has pos- 
itive parity. The hyperfine structure in the excited state 
is the same as in the ground state, but is much smaller 
and is not shown in the main part of the figure. 

The excited state can decay to any of the vibrational 
states of X, but calculations indicate that decay back to 
the v = state is the most probable by far [8j . We use Z 
to denote the Franck-Condon factor for this decay route. 
Within this channel, the decay may proceed either to the 
N = 0, J=1/2,F=0,1 states, labelled (a) and (b) in the 
figure and known as the Qn(l/2) transition, or to the 
N=2, J = 3/2, F=2, 1 states, labelled (c) and (d) in the 
figure and known as the Pi2(3/2) transition. Because the 
hyperfine interaction is not diagonal in J, the two F = 2 
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FIG. 3: Relevant energy levels of CaF. Frequency intervals 
are shown in MHz and are taken from [Tol ] and [111 ]. The plane- 
polarized laser drives the A-X(0-0)Qn(l/2) transition. The 
four magnetic-hyperfine components of the laser transition 
are shown in detail in the inset. The excited state can decay 
back to the v = ground state on any of the branches (a)-(e), 
with branching ratios given by Eq. (IA10|) . Decay to other 
vibrational states of X are also possible and are indicated by 

> o). 



levels are of mixed J character. It follows that decay 
to the upper F = 2 level is also possible, and we label 
it (e). Note that decays to N = 1 are forbidden since 
these states have negative parity, so we have not shown 
these levels in the figure. The decay to the F = 3 level is 
forbidden by the AF = 0, ±1 selection rule. Transitions 
to all other vibrational states are labelled (v > 0), and 
have branching ratio 1 — Z. The branching ratios for the 
decay channels (a)-(e) are calculated in Appendix |A"1 

Taking the quantization axis along the polarization di- 
rection (z), the laser drives only those transitions shown 
in the inset of Fig. [3J We use the index i to denote any of 
the four sublevels of the ground state, and the index j to 
denote the corresponding sublevels in the excited state. 
The coherences reach a steady-state on the timescale of 
the excited-state lifetime. Since the interaction time is 
far longer than this, we can neglect the initial transients 
in the coherences. In this limit, the problem reduces to a 
set of rate equations with the laser excitation rate given 
by: 
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the saturation intensity and is defined by the relation 
I /I s = 20 2 /r 2 . It is the characteristic intensity for the 
first saturation mechanism discussed in Sec. HI 

The Rabi frequency for the transition is Q = EoZij/h, 
where z^ — (i\d.z\j), the laser field is E — EqZ cos(wl£), 
d is the electric dipole operator and £ is a unit vector 
along the z-axis. Using the relations I = l/2eocEo 2 , and 
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we obtain 
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Here, the index k denotes all possible final states in the 
decay of the excited state. We have used the fact that 
the radiation is isotropic (App.[XJ and we have made the 
approximation that u)jk — 2ttc/X, for all final states k. 
This is a good approximation for all relevant states within 
v = whose transition frequencies only differ by a part in 
10 4 . It is less good for the v > decays whose frequencies 
differ by several percent, but since the values of Zjk for 
these transitions is known to be very small @ , the error 
contributed to T by this approximation is negligible. We 
divide the states k into those within the ground state 
J =1/2 manifold, labelled by i, and all other states. With 
this notation, we can write ^X)fe z J fe2 = /Ci^j* 2 ' which 
defines the branching ratio, r, to the J = 1/2 manifold. 
We now use the fact that, for each excited state sublevel 
labelled by j, there is only one non-zero matrix element 
in the sum over i, as illustrated in the inset of Fig.[3J We 
thus obtain a convenient expression for the saturation 
intensity, 
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where we have taken the branching ratio r = 2Z/3 from 
Eq. ifATO]) . 

We now write down the rate equations that govern how 
the populations in the various states shown in the inset of 
Fig. [3] evolve with time. The laser intensity is sufficiently 
low, and its frequency sufficiently far detuned from all 
the other states, that their populations do not need to be 
included in the rate equations. There is a symmetry be- 
tween the (1,-1) and (1,1) states so we do not need sep- 
arate rate equations for each. We denote the populations 
in the (0, 0) and (1, 0) ground states by A^o and N g i, the 
sum of the populations in the (1,1) and (1, —1) ground 
states by N g2 , and the ground state hyperfine interval by 



Here, T is the spontaneous decay rate of the upper level, A g . The notation for the excited state populations and 



S = ujl —^ij is the detuning of the laser angular frequency 
uil from the resonance angular frequency toy, f2 is the 
Rabi frequency, / is the laser intensity, and I s is called 



hyperfine interval is analogous. With the laser detuning, 
S, referenced to the (F' = 1, M' F = 0) <— (F = 0, M F = 0) 
transition frequency, the coupled rate equations for a 
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molecule of speed v are 
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v dNgi_ = R ( S+ A g + A e )[N e0 -N sl ] + ^T[N e0 +N e2 \, 



dz 
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Normalizing all populations to the total population ini- 
tially in the J = 1/2 ground state manifold, the initial 
conditions at the source are N s o = N g i 

N c0 = N cl = N c2 =0. 

Our model does not include the effect of the back- 
ground magnetic field in the experiment. Using a flux- 
gate magnetometer this field was measured to be B x = 
-45±5mG, B y = -212±2mG and B z = -143±lmG. 
The resulting Zeeman splitting between the Mp = ±1 
sublevels is smaller than 1 MHz in both ground and ex- 
cited states. In the ground state, this is too small relative 
to the hyperfine interval for there to be any significant 
mixing of the two hyperfine components. It is also too 
small relative to the natural linewidth to have any in- 
fluence on the excitation rates. The spin precession of 
the excited F' = 1 state is too slow compared to the ex- 
cited state lifetime to have any influence. The excited 
state hyperfine interval, A e , is not yet determined. If it 
is comparable to T, it is much larger than the Zeeman 
splitting, and so the mixing of the F = l, states is small; 
if it is much smaller than T the time evolution of the ex- 
cited state is again too slow to have any influence. These 
considerations show that our neglect of the magnetic field 
is justified. 

The solution of the rate equations depends on the three 
molecular parameters to be determined, the decay rate T, 
the Franck-Condon factor Z, and the excited state hyper- 
fine interval A e . Note from Eq. ((4]) that I s is determined 
once r and Z are known, and so does not constitute an 
extra unknown parameter. To solve the rate equations, 
we also need to know the z-dependence of R, which enters 
via the intensity profile of the laser. This was carefully 
monitored during the experiments. For the model, we 
describe the laser profile by the function 



y) = h 



ta-ah(r/d+L/d) - taah(r / d— L/ d) 
2tanh(L/d) : 
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where r = \J z 2 + e 2 y 2 , L and e describe the size and 
roundness of the beam, and d describes the slope of the 
profile at the edges. This function looks like a top hat 
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FIG. 4: Ground and excited state populations for a molecule 
traveling along the beam axis at 600 m/s and passing through 
a laser beam resonant with the transition from F = 0. Solid 
lines show the (0,0) populations, dashed lines the (1,0) pop- 
ulations, and dotted lines the (1,±1) populations. The grey 
area indicates the profile of the laser beam. The parameters 
of the laser are L = 1.5 mm, d = 0.4 mm and e = 1, and the 
molecular parameters are chosen as Y = 2n x 8 MHz, Z — 0.96 
and A e = 0. (a) and (b) show, respectively, the ground and 
excited state populations when Jo = 0.02/ s . In (c) and (d) 
they are shown for Jo = 4J S . 



when d -C L and like a Gaussian when d ~ L, and we use 
it to fit all of our measured beam profiles. They typically 
have L ~ 1.5 mm, d ~ 0.4 mm and e ~ 1. 

Figure [4] shows how the populations change with z as a 
molecule passes through the laser beam centred at z = 
and resonant with the transition from F = 0. In (a) and 
(b) the laser is weak, Iq — 0.02/ s . The ground state (0, 0) 
population falls as the molecule passes through the beam 
(solid line in (a)), and some small population appears in 
the (1,0) component of the excited state (dashed line in 
(b)). Some of this excited state population decays to 
the (1,±1) ground states, and so their population rises 
(dotted line in (a)). The laser is too weak to excite the 
F = 1 states, whose transition frequencies are detuned by 
many natural linewidths. In (c) and (d) the laser is much 
stronger, Jo — 4J S , and so the population in the resonant 
(0, 0) state falls to zero as soon as the molecule starts 
to enter the beam. Again, this excitation is reflected 
in the population of the (1,0) excited state (dashed line 
in (d)), and some of this decays to the (1,±1) ground 
states as before, resulting in an initial increase in their 
populations. In addition, the laser is now intense enough 
that it can excite population out of the F = 1 levels, 
even though they are far detuned. Consequently, their 
populations also fall, though at a much slower rate than 
that of the F = level. This off-resonant excitation is 
reflected in the excited state populations shown in (d) 
which all become non-zero and follow the profile of the 
laser. The excited state population is greatest in (1, ±1) 
(dotted line), simply because the corresponding ground 
states have more population. The small pedestal in the 
(1, 0) excited state population occurs because population 
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excited first to (1, ±1) can decay to the (0, 0) ground state 
and then be rapidly re-excited to (1,0) by the resonant 
laser. 

The number of photons scattered by a molecule is 
found by integrating the excited state populations over 
time, and multiplying by T. As shown in App. [XJ pho- 
tons are emitted isotropically in the spontaneous decay 
of all four excited states so the signal at the detector 
is directly proportional to the scattered photon number 
without any further complication. In the example dis- 
cussed above, the average number of photons scattered 
per ground-state molecule is 0.28 when I = 0.02J S and 
0.82 when / = AI S . By contrast, when the laser is instead 
tuned to the F = 1 resonance, these values change to 1.07 
and 1.89 respectively. To verify the accuracy of our rate 
equation approach and the neglect of the magnetic field, 
we reproduced the results shown in Fig. 0] using a full 
density matrix calculation that includes all the relevant 
coherences in the problem. The results obtained by the 
two methods were in excellent agreement. The full cal- 
culation reduces the scattered photon number by 0.4% 
in the low intensity case, and increases it by 0.1% in the 
high intensity case. These differences have a negligible 
effect on our final results. 

To model the experiment fully, we first account for 
the divergence of the molecular beam in both the hor- 
izontal and vertical planes. A molecule that arrives at 
the laser at position (x. y) will see the laser frequency 
Doppler-shifted by SiOdopp = ~2irvx/(\D), and will see 
less power than a molecule at y — because of the laser 
profile. Molecules with |x| > 2.25 mm do not contribute 
to the signal since they are outside the detection region, 
while molecules displaced too far in the vertical direction 
see hardly any power and so also contribute negligibly to 
the signal. We solve the rate equations to determine the 
fluorescence signal for many different values of x and y 
within these boundaries, and average the results to ob- 
tain the final signal at the detector. Repetition of this 
procedure for many different values of 8 results in a model 
spectrum. To this we fit a double Lorentzian so that 
it is analyzed in exactly the same manner as the corre- 
sponding experimental spectrum. A set of model spectra 
obtained for a range of laser intensities constitutes a com- 
plete simulation of the experiment. Comparison of the 
experimental data to the simulation data for various T, 
Z and A e , determines these three quantities. 



IV. ANALYSIS 

Figure [5] shows the measured peak heights of the F = l 
and F — spectral lines as a function of laser intensity, 
together with the corresponding results of the model. As 
the intensity is increased from zero, there is initially a 
very rapid increase in the amplitudes of both lines. In 
this regime, the excitation rate, R, is low, so that RT < 1 , 
where T is the transit time of molecules through the 
probe beam. At higher laser intensities, the rapid in- 
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FIG. 5: (Color online) Peak heights of F = and F = l lines 
determined by fitting two Lorentzians, as a function of Jo, 
the laser intensity at the centre of the probe. The other pa- 
rameters of the probe beam are d = 0.38 mm, L = 1.18 mm, 
e = 1. Filled circles: experimental data for F = 0. Open 
circles: experimental data for F = l. Solid lines: correspond- 
ing results of the model, using the best-fit values Z = 0.987, 
r/(2?r) = 8.29 MHz, A e = 4.8 MHz. The inset shows the low 
intensity region in greater detail. Data from several experi- 
mental runs have been combined. 



crease of the amplitudes ceases because all the resonant 
molecules have decayed to other states not resonant with 
the laser. As seen most clearly in the inset, where only 
the low intensity region is plotted, this saturation occurs 
at lower intensity for the F — molecules because, af- 
ter excitation, they are less likely to return to a resonant 
state than is the case for F=l. Inspection of Eq. (|A10[) 
shows that the probabilities are 2Z/9 for (0, 0), 4Z/9 for 
(1,±1) and 6Z/9 for (1,0). For this same reason, the 
amplitude ratio reaches a maximum value of about 4.5, 
as already discussed in the context of Fig. O It is also 
clear from Fig. [S] that the peak heights have not com- 
pletely saturated at high intensity, but instead continue 
to increase with laser intensity, though with a far smaller 
gradient. This is due to the off-resonant excitation of the 
other hyperfinc component, for which the excitation rate 
is small over the entire intensity range shown in the plot. 
This has a greater influence on the F = line, there being 
more F = l molecules to excite off-resonance, so the frac- 
tional rate of increase is larger for this line. At very high 
intensities, approximately 20 times higher than the max- 
imum intensity shown in the plot, the amplitudes stop 
increasing and complete saturation has occurred. At this 
point, molecules in both hyperfine components scatter 
photons until they decay elsewhere, irrespective of which 
hyperfine component is resonant with the probe. Then, 
the mean number of photons scattered per ground state 
molecule would be 1/(1 — r) = 3 if Z were exactly 1. 

Figure [5] shows that the model reproduces the ampli- 
tude data very well over the entire range of intensities. In 
fitting the model to the data, the values of Z, T and A e 
are free parameters to be determined. We explain later 
how we determine these. There is also one additional 
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FIG. 6: (Color online) Widths (FWHM), determined from 
two-Lorentzian fits, of (a) F = and (b) F = 1 lines, as a 
function of 7o, the laser intensity at the centre of the probe. 
The other parameters of the probe beam are d — 0.38 mm, 
L = 1.18 mm, e = 1. The solid lines show the results of 
the model, using the best-fit values Z — 0.987, T/(2n) — 
8.29 MHz, A e = 4.8 MHz. 



free parameter, an overall scaling of the peak heights, 
which is not determined in the experiment since we do 
not measure the absolute flux of molecules. Both the 
F = and F — 1 peak heights are scaled by the same 
value so as to find the best fit to the model. No scaling 
has been applied to the intensity axis. It is worth noting 
that all the details of the model we have outlined need 
to be included in order to obtain good agreement with 
the experimental data. For example, if all molecules are 
assumed to see the same laser intensity, rather than aver- 
aging over the known intensity profile, the 'corner' where 
the on-resonance excitation saturates is far sharper in the 
model than we observe. It is also essential to generate 
complete synthetic spectra, and then fit to these, rather 
than simply determining the amplitudes at the two res- 
onance frequencies. 

Figure [S] shows how the widths of the two lines depend 
on the probe intensity. The observed power broadening is 
related to the saturation of the laser induced fluorescence 
discussed above: as the intensity is increased, the wings 
of the lines increase in amplitude more rapidly than the 
centres, and the lines broaden. It is immediately evident 
in the figure that the F = widths are larger than those of 
F = 1, and that this difference increases with increasing 
intensity. This is because the onset of saturation occurs 



at lower intensity for F — as discussed above. The figure 
shows that the model again agrees very well with the 
data for the particular values of Z, T and A e used here. 
There are no additional free parameters; in particular 
no scaling is applied to either ordinate or abscissa. We 
have only shown the data obtained from one experimental 
run in the figure. Importantly, the widths obtained by 
scanning the AOM frequency (not shown) are consistent 
with the data shown here, indicating that our procedure 
for linearizing the scans is accurate. 

We turn now to the evaluation of the molecular param- 
eters. The dependence of the peak heights on intensity is 
primarily sensitive to the value of Z, as we would expect 
since this determines the mean number of photons scat- 
tered per molecule. Sensitivity of the amplitudes to T 
arises via I s (see Eq. (J5J), and through the degree of off- 
resonance excitation at high intensity. However, given 
the freedom to scale the amplitudes, the sensitivity to 
r is rather weak. The widths, on the other hand, are 
determined primarily by T and are only weakly sensi- 
tive to Z. The excited state hyperfine interval, A e , is 
unresolved in all our spectra, but it does broaden the 
F = 1 line while leaving the F — linewidth unchanged 
(see the inset to Fig. [3|). Furthermore, a larger value of 
A e slightly decreases the gradient of the F — 1 satura- 
tion curve at low intensity, while leaving the correspond- 
ing F = curve unaltered. This happens because the 
laser ceases to be exactly resonant with all three sub- 
components of F = 1 simultaneously. Thus both the 
width data and the peak height data have a small de- 
gree of sensitivity to A e . To determine the central values 
of all three parameters, we simultaneously fit the ampli- 
tudes and widths of both lines to the model, weighting 
each data point by the inverse square of its standard er- 
ror, and so find the parameters that minimize \ 2 m this 
simultaneous fit. This procedure gives Z = 0.987 ±0.012, 
r/(27r) = 8.29 ±0.08 MHz and A e = 4.8 ± 0.4 MHz, where 
we have, so far, quoted only the statistical uncertainties. 
In addition to this simultaneous fit, we have also analyzed 
the amplitudes and widths separately. Fitting to just the 
amplitudes gives Z = 0.98 and lacks sensitivity to the 
other two parameters as discussed above. Fitting to just 
the F' = widths gives T/(2ir) = 8.0MHz, whereas the 
fit to the F' = l widths gives r/(27r) = 8.5 MHz. 

To complete the analysis, we consider the systematic 
uncertainties. There is a 5% uncertainty in the cali- 
bration of the power meter used to measure the laser 
power. Systematically shifting the intensity values, we 
find that this translates to ±0.2 MHz in r/(27r) but does 
not significantly alter the other two parameters. The 
zero-offset in the power measurements, and its uncer- 
tainty, are negligible. Drifts in the molecular flux can re- 
sult in a systematic error in the ratio of the two resonant 
peak heights. We used the frequency-modulated data to 
estimate the size of the systematic drift of beam intensity 
between measurements of the two hyperfine components, 
and find it to be consistent with zero and less than 1.5%. 
This upper bound to the drift contributes uncertainties 
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of _q'oi5 t° ±lMHz to A e and a negligible amount 
to r. Another important uncertainty is the contribu- 
tion of Doppler broadening to the measured linewidths 
arising from the divergence of the molecular beam. The 
Doppler broadening is accounted for in our determination 
of r since our model integrates over the molecular beam 
divergence. In the experiment, an aperture in front of 
the photocathode defines the visible portion of the beam, 
and hence the beam divergence. An aperture placed in 
the molecular beam would have been better. We rely on 
the imaging property of the detection optics to define an 
effective aperture of 4.5 mm in the x-direction. To this, 
we assign a 20% uncertainty. Modelling this, we find 
that decreasing the effective aperture by 20% slightly de- 
creases the linewidths over the entire intensity range, but 
since the Doppler broadening is smaller than the natu- 
ral linewidth the best-fit value of T / (27r) is only reduced 
by 0.2 MHz. Increasing the effective aperture by 20% in- 
creases the linewidths at low intensity as expected, but 
decreases them at high intensity. This is because, on av- 
erage, more power is needed to reach a given level of satu- 
ration when the divergence increases, and so the degree of 
power broadening is reduced. The overall effect is again 
to reduce the best-fit value of T/{2ir) by 0.2 MHz. A 10% 
systematic uncertainty in the diameter of the probe beam 
leads to ±0.14 MHz uncertainty in r/(27r) and a negligi- 
ble error in the other two parameters. Errors introduced 
in linearizing and calibrating the frequency axis can be 
divided into random errors, which are already included in 
the error bars on the individual data points, and a sys- 
tematic error. Based on the consistency with the data 
obtained by scanning the AOM frequency and the ac- 
curacy we obtain for the known ground-state hyperfine 
interval, we estimate that the latter is negligible. 

Combining the error estimates, we obtain the following 
final results: Z = 0.987 ^019' r/(27r) = 8.29 ± 0.30 MHz 
and A e =4.8±1.1 MHz. From T, we find the excited state 
lifetime to be r = 19.2 ±0. 7ns. Our results are consistent 
with the theoretical values, Z = 0.964 and r = 19.48 ns, 
given in [8j . Our measurement of the lifetime also agrees 
with the previous measurement Q, 21.9 ± 4 ns, and is 6 
times more precise. 



V. CONCLUSIONS 



profile and the Doppler shifts of off-axis molecules. Us- 
ing our central values for Z and T in Eq. ([3]) we find the 
saturation intensity relevant to the first saturation mech- 
anism, I s — 222 W/m 2 . Referring back to Fig. [5] we see 
that (apart from the off-resonant excitation) the onset 
of saturation occurs at much lower intensity than this, 
confirming that it is indeed the second mechanism that 
is primarily responsible for saturating the fluorescence. 

Our model can be used as the starting point for mod- 
elling more elaborate experiments that use multiple laser 
frequencies to access several states simultaneously, in or- 
der to scatter more photons. In this context it is in- 
teresting to speculate on the prospects of laser cooling a 
molecule with a structure like that of CaF. Figurc[3]shows 
that 5 laser frequencies (transitions (a)-(e)) are required 
to close off all the transitions within the v = man- 
ifold. For CaF, this would require two separate lasers 
to bridge the 62 GHz gap between N = and N = 2, 
and three AOMs to bridge the hyperfine intervals. Note 
that the system may pump into a dark state unless this 
is destabilized by applying a suitable magnetic field or 
modulating the laser polarization [12| . Note also that, 
as more 'leaks' are shut off, the laser intensity required 
to saturate the transition increases in proportion to the 
scattered photon number, until the intensity reaches I s 
at which point the first saturation mechanism takes over. 
Only under these conditions will the molecules scatter 
the maximum possible number of photons, which for our 
central value of Z = 0.987 is 77 photons. If this could 
be achieved, each molecule could be detected with near 
unit efficiency, since the combined efficiency of PMT and 
collection optics is typically 2%. This photon number is 
still very far from the number required for efficient laser 
cooling - for that, a value of Z much closer to unity is 
required, or further laser frequencies are needed to close 
off the remaining leaks to v > 6]. Speculating a lit- 
tle further, we note that polar molecules such as CaF 
can be brought close to rest using a Stark decelerator 
13L [Tl ] and then confined in a trap for several seconds 
15j. In that case, molecules will already be at temper- 
atures in the 10-100 mK range and only a few hundred 
photon scattering events are needed to reach ultracold 
temperatures. 



By studying in detail the power broadening and sat- 
uration of fluorescence on the A-X(O-O) Qn(l/2) transi- 
tion in CaF, we have measured the excited state lifetime, 
the Franck-Condon factor, and the hyperfine interval in 
the excited state. Our results confirm the accuracy of 
the theoretical work of [8|. Our technique is rather gen- 
eral and can be used to measure the same parameters 
in other molecules. To determine the molecular param- 
eters accurately, the data needed to be compared with 
synthetic spectra generated using a detailed rate model 
of the experiment. The model needs to include all the 
relevant states, and to take into account the laser beam 
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APPENDIX A: BRANCHING RATIOS AND 
ANGULAR DISTRIBUTION OF SPONTANEOUS 
EMISSION 



This can be further reduced by expressing the operator 
in the molecule-fixed coordinate system following [T(| , in 
particular section 5.5.5, to reach 



We calculate the branching ratios for the various decay 
channels shown in figure [3] The calculation is done in 
a basis of states labelled |A, S, S, f2, J, I, F, Mp), usually 
called the (ap) basis and defined in Sec. lIIIl in the context 
of the A state. The excited state is the one of negative 
parity. Omitting the quantum numbers beyond J, this 
state is 



i 

V2 



-41-1 



In X, the expansion of the two states of interest in the 
case (a@) basis is found from Eq. (6.149) of [l6| : 



I AT = T=±) - -J-IO - - - -) + —\0 - —- —- -) 

U ' J 2/ ' ' 2 ' 2 ' 2> 2/ ~ ' ' 2 ' 2' 2 > 2 / > 

I AT = 2 7 = S\ = J_|o 1 1 H\ i In 1 —I —I 2\ 

JV 'i J 2/ \/2 I ' 2 ' 2 ' 2 ' 2 / V2 1 '2' 2' 2'2'' 

Next we deal with the coupling of the electronic an- 
gular momentum and the nuclear spin. All the relevant 
states are diagonal in J apart from the two states that 
have F = 2, as discussed in Sec. IIIII Diagonalization of 
this part of the problem can be done exclusively within 
the N = 2 manifold following the procedure given in the 
appendix of [13] and using the spin-rotation and hyper- 
fine parameters given in [1 01 ] - Labelling the upper and 
lower F = 2 levels by |2±), the result is 1 



|2+> = ||J= §,/=!,*■= 2}+ &\J=ll=±,F = 2), 
|2_) = ^|J=|,J= \,F = 2)-\\J=\j=\ i F = 2). 

Now we calculate the matrix elements of the electric 
dipole operator, written in spherical tensor notation as 
T^(d), between basis states of A and X. Applying 
the Wigner-Eckart theorem and noting that the opera- 
tor acts on the electronic angular momentum, the matrix 
elements factorize to 



m fi = (f\ T p 1 \d')\i) = mi m 2 m 3 
where p = 0, ±1 labels the polarization, and 



(Al) 



mi = (-1) 



F'-M' 



(A2) 



F' 1 F 
-Mp p M F 

(-l)^'+'+V(2F' + l)(2F+l)j 

(A3) 

™ 3 = (A', S, £', J'| (d) 1 1 A, S, S, O, J). (A4) 



m 2 



1 That the value of the coefficient, numerically 0.4995, should be 
so close to i seems to be coincidental. 



= E(- i ) j '- n V(2J'+i)(2j + D( J QI \ i 



A',5,E'|T g W(d)|A,S,E\ (A5) 



m 3 



To obtain the branching ratios, we do not need to 
evaluate the final factor in the above equation since it 
is common to all the branches. Note that the electric 
dipole operator cannot change either the spin or its pro- 
jection, so the only non-zero matrix elements in the above 
sum are those that have S = E'. Since A = whereas 
A' = ±1, it follows that ft = Q' ± 1 and so the q = 
term in the sum is zero. Writing the transition matrix 
elements in the form m^fl' , fi), and inserting the val- 
ues for all the quantum numbers fixed in the problem 
(|A'| = 1, A = 0, |E'| =S=LJ' = I), we have 



m 3 (L |) = to 3 (-3, -\) = 0, 



\/2(2J + 1) 



4 i 



i j 
i 

2 



TO4 



m 3 (-I,i) = (-l)2-V3(i,-i) 
m 4 = (l,|,-i T +1 (1 \d) 



(A6) 
(A7) 

(A8) 
(A9) 



Putting all this together, we get the branching ratios 
for the entire set of decay channels. For reference, the 



'ical results 
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(A10) 



Here, columns are labelled by the (F 1 , Mp) values in the 
upper state, and rows labelled by the (N, F, Mp) values 
in the lower state. 
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The angular distribution of the spontaneously emitted 
photons is proportional to sin 2 9 for AMp — {) transitions 
and to (1 +cos 2 #)/2 for AMp = ±1 transitions, where 
9 measures the angle between the observation direction 
and the z-axis. For the decay of the (1,0) and (0,0) 
states, Eq. (|A10p shows that the total branching ratios 



for AMp — 0, ±1 are each 1/3 and so the emission is 
isotropic for these states. For the (1, 1) state the branch- 
ing ratios are 1/2,1/3,1/6 for AMp = -1, 0, +1 transi- 
tions and so the emission from this state is also isotropic. 
For the (1,-1) state these branching ratios are reversed 
and the emission is again isotropic. 
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